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I. INTRODUCTION 



and the progresses in 



The study of neutrino mixing in the context of quantum field theory (QFT) 
the understanding of such a phenomenon 0, Q| , together with the definitive experimental proof of neutrino oscillations 
3) [3 1 opsn new scenarios to the research in fundamental physics. Indeed, it has emerged an unexpectedly rich non- 
perturbative structure associated to the mixing of neutrino (as well as boson [ll]) fields, hidden in the vacuum for the 
flavor fields. This has been shown to be a condensate of massive neutrino-antineutrino pairs. Several consequences 
of this discovery have been analyzed, including the effects on flavor oscillation formulas [l|, |5| and the implications in 
cosmology and astro-particle physics 12, llSl, \\A \ . 

In this review we show that the energy content of the neutrino mixing vacuum condensate 
component of dark energy 13| that, at present epoch, assumes the behavior and the value of the observed cosmological 



13 | can represent a 



constant. We compute such a value and show that, above a threshold, it is slowly diverging and its derivative with 
respect to the cut-off value goes actually to zero, which allows to use the cut-off at its electroweak scale value, provided 
one limits himself to considering the two (lighter) generations in neutrino mixing. In the case one includes mixing 
with heavier neutrinos, a value of the dark energy compatible with its upper bound is obtained for a cut-off of the 
order of the natural scale of neutrino mixing. As we will show in a forthcoming paper [15|] , the use of such an infrared 
momentum cut-off is motivated, at the present epoch, by the negligible breaking of Lorentz invariance due to the 
vacuum condensate caused by neutrino mixing . 

The remarkable improvement of the dark energy value computed in the present paper, with respect to the disagree- 
ment of 123 orders of magnitude in standard approaches [l^ . makes the present treatment worth to be discussed. 



17| . The link comes from the neutrino- 



Our result links together dark energy with the sub-eV neutrino mass scale 
antineutrino pair vacuum condensate. 

We point out that our work differs from the approach [ist based on vacuum contributions from light particles 
like neutrinos and axions. In the present review we do not resort to axion contributions. Nevertheless, it is in our 
future plan to compare the two approaches in order to clarify the differences and the similarities. Moreover, the 
non-perturbative feature here presented leads us to believe that a neutrino-antineutrino asymmetry, if any, related 
with lepton number violation 3| would not affect much our result. Also this point deserves to be better clarified. 
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which we plan to do in a subsequent work. Finally, we observe that the non-perturbative contribution discussed in 
the present review is of different origin with respect to the vacuum energy contribution of massive spinor fields arising 
from a radiative correction at some perturbative order [2^. 

The review is organized as follows. In Section II, we introduce the neutrino mixing formalism in QFT. In Section 
III we present the neutrino mixing contribution to the dark energy of the universe. Conclusions arc drawn in Section 
IV. 



II. NEUTRINO MIXING IN QUANTUM FIELD THEORY 



The neutrino mixing phenomenon was firstly studied in the context of quantum mechanics (QM) j2lL l22l . l23l . \24 



25 



26 



27 



and subsequently analyzed in the framework of the QFT formalism I, which we 

shortly summarize in the following (for a detailed review see 

In the very effective Pontecorvo's formalism neutrino mixing is considered from the standpoint of QM and the 
attention is focused on the mixing of " states" . The fact that neutrinos are actually described by field operators is 
completely neglected. The reason for that is the necessity of the effectiveness of the formalism which is required 
to readily fit the experimental search for mixing and oscillations. Therefore, any simplification of the matter to be 
treated is adopted, provided the resulting description would be sufficiently accurate and descriptive/predictive of 
the experimental observations. As a matter of fact, the successive development of the experimental search has been 
supporting such an attitude. From a theoretical point of view, there is, however, the necessity to understand how 
mixing and oscillations can be properly described in the realm of QFT, which provides anyway the proper setting for 
neutrino dynamics, as known since the birth of QFT. It is then also necessary to understand how the correct formalism 
connects to the Pontecorvo's approximate scheme. This has been indeed the program of the research line which has led 
to the QFT formulation of the neutrino mixing and oscillation. Such a program has been successively extended so to 
incorporate other particle mixing (quark mixing, boson mixing) . Here we only sketch the skeleton of the QFT mixing 
formalism and to do that we consider two neutrinos. Extension to three neutrino (to any number of generations, in 
principle) is in the literature [5]. The reader who wants the guaranties offered by a rigorous mathematical proof of 
our treatment may usefully read the papers in Refs.^j. 

The mixing transformations for two Dirac neutrino fields are 



Veix) — i^iix) COS0 + 2^2(2;) sin 9 
v^ix) — —vi(x)sniO + V2{x)cost 



(1) 



where Ve{x) and v^{x) are the fields with definite fiavors, 9 is the mixing angle and Vi and V2 are the fields with 
definite masses mi 7^ m2'. 



1,2, 



(2) 



k,i' 



1,2, r = l,2 



with a^^S) = "k.» 6^^'^'=.'*, /3k^,(t) = I3\^^^ e*"'=.'*, and ujk,t = + mf. The operators aj; and /3j 
annihilate the vacuum state |0)i.2 = |0)i ® |0)2: akJ0)i2 = /3J^j|0)i2 = 0. The anticommutation relations are: 
{vf{x),v'^j\v)^ = ^^(x-y)(5„/3(5,y, with a,(3 = 1,...4, and |ak_,,aqj j = &]^q_&rs5if, {^k,*'/^qj} = &\^q.5rs&i3, 



with i,j = 1,2. All other anticommutators are zero. The orthonormahty and completeness relations are: u^iU^^ = 

The mixing transformation Eqs.([T]) can be written as [ij: 

,y:{x)^G,\t),.n^)Ge{t) (3) 
'^:{x)=G,\t),.^{x) Ge{t) 



where the mixing generator Gg{t) is given by 

Gg{t) — exp 



d^X (l'l{x)l'2{x) — l'2{x)'^l{x) 



(4) 



At finite volume, Ge{t) is an unitary operator, Gg^{t) = G-g{t) = Gg{t), preserving the canonical anticommutation 
relations; Gg^{t) maps the Hilbert spaces for i^i and V2 fields Tii_2 to the Hilbert spaces for flavored flelds Tie,ii'- 
Gg^{t) : 7ii_2 'He.fA.- In particular, for the vacuum |0)i^2 we have, at finite volume V: 

\m)e,^^Gg\t)\Q),,2. (5) 

|0)e./j is the vacuum for 7ie,p, which we will refer to as the flavor vacuum. The explicit expression for \Q)e,ti at time 
t = in the reference frame for which k = (0, 0, |k|) is 



r,k 

+ e'-sin^ e |^k||C/k|(<i/3:li - «L^2/5:l2) + sin' ^ l^^kpaL!i/?1.2«L^2/3l.il |0>i,2 ■ (6) 



Eq.® exhibits the condensate structure of the flavor vacuum |0)e,p. The important point is that i,2(0|0(<))e_^ 0, 
for any t, in the infinite volume limit Thus, in such a limit the Hilbert spaces 7ii,2 and 'He,^^ turn out to be 
unitarily inequivalent spaces. We remark that \Q)e,^^ is the physical vacuum as we will see below. 

In the QM formalism we " cannot" have two unitarily inequivalent Hilbert spaces for the simple reason that the von 
Neumann theorem forbids the existence of unitarily inequivalent representations of the canonical (anti-)commutation 
rules whenever the system has a finite number of degrees of freedom, such as those in QM. It is quite obvious that 
this cannot be the case for neutrinos. Since they are quantum fields, by definition they are described by infinitely 
many degrees of freedom and thus von Neumann theorem does not hold. This point introduces a crucial difference 
between the QFT formalism and the QM approach. 

The fiavor annihilators, relative to the fields Ve{x) and v^ix) at each time, are given by (we use (cr, i) = (e, 1), (/i, 2)): 

^ G'g\t) al^,{t) Gg{t) , 
PlAt) ^ G-\t) Pl,{t) Ge{t) . (7) 

The flavor fields can be expanded in the same bases as Vi. 



k.r 
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The flavor annihilation operators in the reference frame such that k = (0, 0, |k|) are: 

alJt) = cose al.it) + sine [\U^\ al,{t) + \V^\ pX2(t)) 
al^it) = coseal^it) - sine {\U^\ al,{t) - e'' |^k| /3la(i)) 

pL^jt) = cosepL^^,{t) + sine (|[/k| /^Lk.zW - l^kl (9) 

Ck,pW = cos^Ck.aW - sine (|t/k| /S^k.iW + e'' l^k| , 

where e'' = (—1)'^ and 

|[/k| = Wk!j'"kj = ^^-k,,'^-kj ; * 7^ i ; = ejje'' w^^^wllkj, no summation (10) 

with Eij — 0, 1, —1 for i — j, i < j, i > j, respectively. We have: 

|C/kP + |VkP = l. (12) 
The number of condensate neutrinos for each k is given by 

eAOK]^<Me,^ = e,MPKiPU(^) e.^ = S^nH \Vu\' , (13) 

with i = 1,2. Equivalently, i,2(0|ak!.aLJ0>i.2 = iA0K]aPlJ0)h2 = sin^^ jVkp , with a = e,/i. 

The Bogoliubov coefficient |VkP appearing in Eg. p^ can be written as a fimction of the dimensionless momentum 



and dimensionless parameter follows, 



|2 1 P 



2 



|y(p, a)|^ = - 1 - " ^' ^ ^ . (14) 

2 \ +1) +«P / 

From Fig.l we see that the effect is maximal when p = 1, i.e. for |kp = mim2, the natural scale of the neutrino 
mixing. |VkP goes to zero for large momenta, i.e. for |kp 3> mim2, as |VkP ~ ^'^^i ■ I* acts as a "form factor" in 
the k space controlling the neutrino vacuum condensate. We thus find that the Pontecorvo's formalism is nothing 
but the relativistic limit of the QFT formalism: |Vkp and |?7kP 1 for |kp 3> mim2. In the Pontecorvo's 
formalism the non-perturbative contributions from the vacuum condensate are thus missing. This is the meaning of 
the approximation made in the QM treatment of the mixing. Missing these condensate contributions is of course 
of no relevance for the experimental observation of neutrino oscillations at today instrumentation resolution. These 
contributions might play, however, a relevant role in the study of the cosmological background. The fact that |Vkp 
contributes maximally for low energies suggests indeed to us that the contribution of the mixing phenomenon may be 
taken as a good candidate in the study of dark energy. 

We point out that, since |0)i^2 and |0)e,/x are unitary inequivalent states in the infinite volume limit, two different 
normal orderings must be defined, respectively with respect to the vacuum |0)i^2 for fields with definite masses, 
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Figure 1: The fermion condensation density a)|^ as a function of p for a = 0.98 (solid line) and a = 0.5 (dashed hne). 



denoted as usual by : ... :, and with respect to the vacuum for fields with definite flavor |0)e,/i, denoted by 
The Hamiltonian normal ordered with respect to the vacuum |0)i_2 is given by 

■H: = H-1,2 (0|i7|0)i,2 = H+ 2 [ cP\<i{uka+uJk,i)^Y.Y. I ^'kc^M[<,aM + /^k^J , 

i r 

and the Hamiltonian normal ordered with respect to the vacuum |0)e,^i is 

::H : = H - ,,^(0(t)|ff |0(i))e^^ = H + 2^ d^k{L0k,2 + tOk,i) (l - 2\V^f sin^ 9) . 



(15) 



(16) 



Note that the difference of energy between |0)e,^ and |0)i^2 represents the energy of the condensed neutrinos given in 
Eq.idni) 



(O(i)l : |0(t))e,^ = eAm\H\m)e,^ 



1,2 



m\o) 



1,2 



Asin^d J d3k(c^,,2+^fe,i)|Vu|' , 



(17) 



and gives the "energy gap" protecting the flavored neutrinos from turning into the mixing component neutrinos vi 
and 1/2- In the following we show that the energy of the condensed neutrinos can have cosmological implications, 
indeed it can contribute to the dark energy of the universe. 

Before considering cosmological aspects of neutrino mixing, in order to better understand the meaning of Ea. (|17p . 
let us introduce the operator A{t) that satisfies Eqs. l|Al[ ) - (jA4[) given in Appendix A. By defining the operator 



we have 



H'{t) = :H : -A(t) , 



{KA'^)\H'{t)VKe{'^)) = c^feacos^e +c^,,2sin2 0, 

{^lMH'{t)\yl^{t)) = C.fc,2COs2 +C.fe,iSin2 0, 



(18) 

(19) 
(20) 
(21) 



(22) 
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Eas. (fT9| - (l2ip coincide with the ones obtained in QM by using the Pontecorvo states. Moreover the uncertainties in 
the energy H'{t) of the multi-particle states (jASp . (jA6| are zero such as in QM, are zero the uncertainties in the 
energy H of the multi-particle states. H'{t) is explicitly given by 

r 

+ u,, + <,(i)«Le(i) + /3:le(i)/31k,,(0 + /3:l,,(i)/3Le(i)) ] , (23) 

where Wee = ^^fc,! cos^ 9 +u}k,2 sin^ 9, = ujk^i s\t? 9 +i0k,2 cos^ 9, w^e = {^k,2 — sin 6* COS 6*. From Eq. ([23|l we 
have at any time t 

e,^(OW|if'(0|OW)e,^ = 0. (24) 

Thus we obtain 

e,^(0(i)| : H : |0(i))e,M ^ e,p(0(t)|A(i)|0(i))e,^ • (25) 

That is, the operator A{t) is the part of the Hamiltonian : H : that give rise to the neutrino mixing condensate. The 
operator ([23|l can be also written as 

H'{t) =::H::-B{t), (26) 

where Bit) satisfies Eqs. (|A7p - (jA10[) presented in Appendix A. 

We remark that one might also consider the "effective" Hamiltonian approach by incorporating into the QM 
treatment the condensate contributions computed by using the operator A{t) (or B(t)). 



III. NEUTRINO MIXING AND DARK ENERGY 



In this Section we show that the energy density of the neutrino vacuum condensate can represent an evolving 
component of the dark energy. The non-zero value of |Vkp for long wavelengths, namely its behavior at very high 
momenta, together with the negligible breaking of the Lorentz invariance of the vacuum condensate at the present 
time, can be responsible of the very tiny value of the cosmological constant. 

Let us calculate the contribution p™" of the neutrino mixing to the vacuum energy density in the Minkowski metric. 
The energy- momentum tensor density T^u{x) for the fields vi and V2 is 

: ^ : (^Mlt. d. *™(a:)) : (27) 

where '!',„ = {i^i,i^2)'^- The symbol : ... : denotes the normal ordering with respect to |0)i.2. 

In the early universe epochs, when the Lorentz invariance of the vacuum condensate is broken, p™" presents also 
space-time dependent condensate contributions. Then p™" is given by the expectation value of the (0,0) component 
of T^,^{x) in the fiavor vacuum \0{t))^ 

pZc = '?oo ■■ T°"iO) : m)),,^ , (28) 



where : Tqq : = : H := J d?x : Tqq{x) : . Note that Too is time independent. We obtain 



K 

pZc = -sin^e / dfcPK,! +c^fc.2)|14|' 

TT 



(29) 



where the cut-off K has been introduced. Exphcitly 

1 . 



mix 
rvac 



2tt 



ml log 



m2 



ml log 



L ^ + mf j I 



(30) 



The contribution of the neutrino mixing to the vacuum pressure is given by the expectation value of 7jj (where 
no summation on the index j is intended) in |0)e,;^: 



where Tjj — J (fixTjj{x) . Being 



in the case of the isotropy of the momenta: T^^ = T^^ = T^^, we have 



Pv 



■ sm 





1 


1 


/ dkk^ 






la 







(31) 



(32) 



(33) 



Explicitly 



mix 
r^vac 



— sin^ 0{m2 — mi)^K (^rn2^ + mf — mi\J~K^ 



2 



arctan 



{2m\ + mim\) log 



n2\/K^ + ml 

1712 



K - 



ml 



arctan 



V 



11^2 — '"'1 

(2m2 + m\m2) log 



K + ^K'^ + m{ 



mi 



(34) 



The state equation of the vacuum condensate is defined as i/;™" = p"ac /p^ac- plotting ly™*^ as function of the 
momentum cut-off K (Fig. 2) we have that u;™" = 1/3 when the cut-off is chosen to be K ^ mi, m2 and u;™" goes 
to zero for K < yjm\m2- 

The neutrino vacuum condensate assume a different behavior at the present epoch. Indeed, the breaking of the 
Lorentz invariance is now negligible and, at present time, p™^ comes from space-time independent condensate con- 
tributions. Then the energy- momentum density tensor of the vacuum condensate is given approximatively by 



T^. : |0) 



^E^ 



(27r)3 ^'^^ 



Vi(x)v^{x) : |0)e,p = Viiu pT 



(35) 



Since 77^1^ = (img(l, — 1, — 1, — 1) and, in a homogeneous and isotropic universe, the energy-momentum tensor is 
= diag{p ,p ,p ,p), the state equation is then p""^ « — P™^, that is, the neutrino vacuum condensate today has 
a behavior similar to the cosmological constant 

2 



13|. Explicitly, we have 



p™" = - sm- 

TT 



K 



in'^e / dkk^ 



(36) 
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Figure 2: The adiabatic index ui™" as a function of cut-off K. 

We observe that the value of the integral ([55]) is conditioned by the appearance in the integrand of the |VkP factor. 
The integral, and thus would be zero for |Vkp = for any |k|, as it is in the usual quantum mechanical 

Pontecorvo formalism. In the present QFT formalism |Vkp is not identically zero for any |k|; it goes to zero only for 
large momenta, getting its maximum value for |k| — y^mlmj {p — 1, cf. Section II), thus maximally contributing 
in this (infrared) region: we thus see that contributions to dark energy mostly come from long wave-lengths, short 
wave-lengths at most producing local dishomogeneities. Proceeding in our calculation, we obtain 

1 9 



Pa = 2^^^"^ ^"^^ 



+ 2 



: arctan 




{ml + 2ml + miml) log ( ^ — '-^ ) + {mf + 2m\ -f m{m2) log ( ^ — ) \. (37) 

To better understand the meaning of Eg. ([57)1 . we report the behavior of p™'^ for K ^ toi, m-i,'- 



Pa 



^ sin^ 6{m, - mM2\^^= arctan f V^EM) _ < arctan f V^EM) 



^ ml — m\ ^ nil ' mf — m^ ^ ^2 

- (rn^ + 2™? + mimj) log ( — ) + (to? + 2to^ + mim2)log ( — ) |. (38) 

Vto2 / \mij) 

This shows that the integral diverges in K as mf log {2K/mj), with i,j = 1, 2. As shown in Fig. 3 the divergence in K 

is smoothed by the factor mf. By using the electroweak scale cut-off: K = 100 GeV, for neutrino masses of order of 

lO^'^ey and Ato^j ~ 7 x 10"''^ey^ we have p™" w 2.9 x 10~'*^Gey^, which is in agreement with the observed value 

of cosmological constant. From Eq. ([38|l one also sees that ''^"J^^^^ cx -> for large K. 

The result we have obtained is quite instructive also since it tells us that the value of |VkP, for any |k|, contributing 

to the observed value of p™*^ is the one related with such a mass scale (the dependence of |Vkp on the masses is 

shown in Eq. (fTTj) . see also Fig.l). The computation of p™" turns out to be sensible to small variations in the values 

of the neutrino masses and of Ato^, these last ones affecting the value of the multiplicative digits of lO^'^^GeV^'*. 
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K(GeV) 



Figure 3: The neutrino mixing dark energy as a function of cut-off K. 

Above we have derived the contribution p™'^ arising from mixing of the two hghter neutrinos. If mixing involving 
heaviest neutrinos are included, the obtained value for dark energy, for a value of the cut-off of order of the electroweak 
scale, turns out about 4 orders of magnitudes higher than the observed value of dark energy. In such a case, a value 
of the dark energy compatible with its upper bound is obtained for a cut-off of the order of the natural scale of the 
neutrino mixing. Such a small cut off on the momenta is imposed by the negligible Lorentz invariance breaking at 
the present epoch, as we will show in a forthcoming paper There we will present also the explicit computation 
in curved space-time and we will show that the mixing treatment here presented in the flat space-time is a good 
approximation in the present epoch of that in FRW space-time. 

IV. CONCLUSIONS 

In this report we have presented the main features of neutrino mixing in the context of quantum field theory and 
we have shown that neutrino mixing may contribute to the value of the dark energy exactly because of the non- 
perturbative field theory effects. In particular, we have shown that, at the present epoch, the vacuum condensate 
generated by neutrino mixing behaves as the cosmological constant. Its observed value is obtained for a cut-off of the 
order of electroweak scale when the two lighter neutrinos are considered and for a cut-off of the order of the natural 
scale of neutrino mixing in the case one includes mixing with heavier neutrinos. 
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Appendix A: EXPECTATION VALUES OF THE OPERATORS A, AND B 

The operator A{t) satisfies the relations: 

i'^LMmWLeit)) = 2^,,i|yk|'sin2 0; {,ylJt)\A{t)\,.lJt)) ^2o.,M'^i^'0, (Al) 
KeWI^WK^W) = K^(t)|A(OK,W> = K2-'^M)(|C^k|-l)sin0cos0, (A2) 

and similar relations hold for the anti-particle states \T'a-{t)), moreover 

Kep^WI^WKeW) = 2e^uj,s sm'e\U^\\V^\; {K.^.M A{t) Wljt)) ^ ~2e^u;k,2 sin^ ^ |Fk|, (A3) 
K,e,M'^(t) KA*)) = K,,MA{t) WUt)) = {cok,2+iOk,i) iMcl sin^ COS0, (A4) 

where, at time t, the multi-particle flavor states are defined as: 

l^k,ee> W> ^ <e W /^l.e W "L^ W \m)e,, , (A5) 
Wl,-,eit)) ^ <pW/3l,^W<eW|0W>e,,.. (A6) 

The operator B in Eg. ipS)) has the expectation values given below: 

KeWI^WKeW) = -2uj,^2\V^\'8mH; {,yl^{t)\B{t)WlJt)) = -2co,MVu\'sm'd, (A7) 
{KAt)\B{t)Wl^it)) = {iyl^it)\B{t)\iylJt))^ic,k.2-iOk,,){\Uk\~l) sine cose, (A8) 

and similar for |;7cr(t)), moreover 

K,e^MBit)WlAt)) = 26'- sm'e\U^\\V^\; {>^L,,M B{t) Wl^{t)) ^ -2e'-c.M sin^ (A9) 
Kef.M Bit) WLJt)) - B{t) = (a;,,2 +c.,,i) 1141 sin0 cos0. (AlO) 
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